Cold atom dynamics in a quantum optical lattice potential 
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We study a generalized cold atom Bose Hubbard model, where the periodic optical potential is 
formed by a cavity field with quantum properties. On the one hand the common coupling of all 
atoms to the same mode introduces cavity mediated long range atom-atom interactions and on the 
other hand atomic backaction on the field introduces atom-field entanglement. This modifies the 
properties of the associated quantum phase transitions and allows for new correlated atom-field 
states including superposition of different atomic quantum phases. After deriving an approximative 
Hamiltonian including the new long range interaction terms we exhibit central physical phenomena 
at generic configurations of few atoms in few wells. We find strong modifications of population 
fiuctuations and next-nearest neighbor correlations near the phase transition point. 
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Laser fields can nowadays be routinely used to create 
tailored optical potentials for ultracold neutral atoms 1] . 
Loading an atomic BEC into such a periodic standing 
light pattern allows to experimentally implement sys- 
tems, which are well described by a Bose Hubbard Hamil- 
tonian with externally controllable parameters Q, |M Q - 
In some recent spectacular experiments the predicted 
Mott-insulator to superfluid quantum phase transition 
has been observed fw]. As the light fields are normally 
intense and strongly detuned from any atomic transition, 
their properties can be safely approximated by prescribed 
classical fields independent of the atomic state. How- 
ever, this is invalid if they are confined within an optical 
resonator. For a sufficient atom number N and atom- 
field coupling g the fields become dynamical quantities 
depending on the atoms. In addition in a high-Q cavity 
the quantum properties of the field get important and the 
atoms move in quantized potentials ^Jj. Ultimately this 
allows different states of the lattice field (e.g. different 
photon numbers) to be quantum correlated with differ- 
ent quantum phases of the atoms. As a striking example 
the atoms could be in a superposition of a Mott insulator 
and a superfluid state connected with a different cavity 
field amplitudes. Even in the classical field limit of high 
photon numbers all atoms see the same field state and 
thus we get new long range atom-atom interactions. In- 
terestingly the idea of implementing such combination of 
cavity QED and a BEC has experimentally made such 
rapid progress recently, that its success can be expected 
very soon Q . 

In this work we discuss the basic physical properties of 
such a generalized model of ultracold atoms in a periodic 
potential generated by a quantized field mode. In a first 
step we derive an approximate Hamiltonian analogous 
to the Bose Hubbard Hamiltonian including a quantized 
potential. Its basic physical implications are then exhib- 
ited in two cases: (a) a strongly damped cavity, where 
the field dynamics can be adiabatically eliminated, which 
leads to a rescaling of the coupling parameters and new 
long range atom-atom coupling terms and (b) a dynam- 



ical model where the field is approximated by its time 
dependent expectation value derived from a dynamical 
equation containing atomic expectation values. 

Model: Let us consider N two-level atoms with mass m 
and transition frequency LUa strongly interacting with a 
single standing wave cavity mode of frequency ujc i^a- 
The system is coherently driven by a laser field with fre- 
quency LUp « LUc through the cavity mirror with ampli- 
tude T]. Alternatively the atoms are illuminated transver- 
sally to the cavity axes (see Fig. ^ with amplitude C- 
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FIG. I: (color online). Sketch of setup 



Including damping, the dynamics is given by the mas- 
ter equation for the atom field density operator g — 
J-[-ff, g]+Cg, where the Liouvillian C models dissipation. 
For large atom-pump detuning spontaneous emission is 
negligible and cavity loss k will be the dominant dissipa- 
tive process, i.e., £g = k (2aga^ — a)ag — go) a). Here a 
is the annihilation operator for a cavity photon. 

As convenient choice we rewrite the Hamiltonian in 
a second quantized form, where the direct interaction 
between the atoms is modeled by a pseudopotential and 
characterized by the s-wave scattering length as'. 



H = j d^a;^'t(x)i7o*(x) 

+ / d3a;*t (x) (x) ^ (x) ^ (x) . (1) 

2 m J 



2 



Here ^ (x) is the field operator for the atoms and Hq 
is the single-particle Hamiltonian in rotating-wave and 
dipole approximation: 



contain more interesting physics, we will concentrate on 
this limit to be able to proceed analytically. Explicitly 
the coupling matrix elements read: 



Hq = hAaa^<J — hAcO^a — ihg(x) (a'^a — (7 a 



ihh (x) C (ct^ — (7 ) — ihrj (a — a^) 



(2) 



where = Up — uja and Ac = ujp — Wc is the atom- 
pump and cavity-pump detuning, respectively. Along the 
cavity axis (x-direction) the atom-field coupling is set to 
g{x) — gocos{kx), while the transverse laser beam forms 
a standing wave with amplitude h (x) = ho cos{kpy) 
in the y-direction, where we set y = for our ID- 
considerations. In the regime of low saturation j2] (large 
Aa) adiabatic elimination of the excited atomic state 
in ^ then leads to: 



Hq = h cos^ikx) (hUoa^a + Vd) — fiAca^a 

— ihrj (a — a^) + hrjcs cos{kx) (a -|- a^) . (3) 

The important parameter Uq = g^/ Aa here is the op- 
tical lattice depth per photon ||6j and also gives the re- 
fractive index of one atom at an antinode. The term 
containing ryofr = gohoC/Aa represents an effective pump 
through atomic scattering into the mode. Along x the 
cavity field forms an optical lattice potential with period 
A/2 and depth hUoa^a. For the sake of generality we add 
an extra classical potential Vd as well. 

To derive a generalized Bose-Hubbard Hamiltonian we 
expand the Bloch states of a single atom inside the lattice 
using localized Wannier functions |lCl| and rewrite the 
field operators in Eq. in these functions, keeping only 
the lowest vibrational state at each site (lowest band) 
* (x) = J2i biW (x - Xi) to get: 



kd k,l 

+ fir]cs (a -I- at) ^ Jk,ib\bi - ifirj (a - a)) 

k,l 

jY.^ibk(bibk-l)-hAca^a. (4) 



The operators 6^ (bk) correspond to the creation (anni- 
hilation) of an atom at site A; and the on-site interaction 



of two atoms is given by U 



/d3a;|(x)r. As the 



nonlinear part of the nearest-neighbor interaction is typ- 
ically two orders of magnitude smaller than the on-site 
interaction it is neglected as usually. In contrast to the 
classical field case |3| the appearance of the cavity field 
operator does not allow to reassemble all hopping terms 
to a single expression. To be still able to proceed analyti- 
cally we assume a weak dependence of w{x) on the cavity 
photon number. Although the opposite limit might even 



Ekd = J d^xw{x-:>ik) (^-^V^^ w(x-xO, (5a) 
Jk,i = j d'^a; w (x — Xfc) cos^(fc.x)w (x — Xi) , (5b) 
Jk,i = / d'^a; w (x — Xfc) cos(fca:)w (x — X;) . (5c) 



These matrix elements are symmetric, i.e., Ek,i = 
Ei^k,Jk,i = Ji,k and Jk,i = Ji,k and the on-site ele- 
ments Ek^k and J^^k are independent of k. As the next- 
nearest neighbor terms are typically two orders of magni- 
tude smaller than the nearest-neighbor amplitudes they 
are omitted too. Note that in the case of transverse 
pumping two adjacent wells acquire different depths since 
the cos in Q changes sign periodically, which implies 
Jk,k — — Jk+i,k+i- The Hamiltonian Q now reads: 



H = EqN + EI3+ {hUoo^a + V^) [JoN + JB 
+ hrjcs (g + flt) Jq y^(-l)^+^rtfc - hAca^a 



ihrj (a 



k 

-E 



nk [fik - 1) , 



(6) 



where we introduced N — rik — ^I^fc (number 
operator) and B = {b\_^_ibk + h.c^ (jump operator). 

i?o, ^0 and Jo are on-site matrix elements, whereas E 
and J are the site-to-site hopping elements. The Hamil- 
tonian is a central result of this work and gives the 
starting point to discuss the physics below. 

Lets first look at the light field dynamics and write 
down the corresponding Heisenberg equation: 



(a, - [/o (Jo^ + JB 
iilosJoY^i-'^f^^fik. 



a + Tj 



(7) 



We see that the quantum state of the field depends not 
only on the number of atoms N but also on coherences 
via 13. Particularly interesting effects can also be ex- 
pected from the last term describing transverse pumping 
as the corresponding operator has an alternating sign for 
neighboring wells. Hence it vanishes exactly for a Mott 
insulator state, while it gives a nonzero contribution for a 
superfluid state. Nevertheless we will concentrate here on 
the more simple setup pumping acts via a cavity mirror 
and set T^off = below. 

Bad cavity limit: The common interaction of all atoms 
with the same field implies a complicated dynamics. 
Luckily in typical setups the field damping rate k is the 



3 



fastest time scale in the system. This allows to eliminate 
the cavity degrees of freedom by formally solving equa- 
tion for a = ry/[K-i(Ac-?7o(Jo^ + -'^-S))] and insert- 
ing this back to I^J. As TV commutes with B this gives 
no ordering problem. For a fixed atom number N = {N) 
we then expand a to second order in the small tunneling 
matrix element J: 
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iUoJ 
K ~ iA' 



B 



:B' 



(8) 



where = Ac — UqJoN is a rescaled detuning, so that 
we have: 



FIG. 3: (color online), (a) Energy difference AE as function 
of Ac for a single atom in two wells for Uo = —1.2k (solid 
line) and Uo = —0.4k (dashed line). The associated lattice 
depth is shown in (b). The other parameters are rj = 2k, and 
= ~4Er. 



H 



(«2 + Af)' 2^ 

Obviously H now contains cavity induced rescaling of 
the tunnel coupling proportional to B as well as nonlocal 
correlated two atom hopping terms proportional to B^ 
mediating long range interactions. In Fig.|21we show the 
excellent agreement of the field calculated from the ex- 
pansion compared to a full numerical inversion within 
a large parameter range as used below. Note that the 
matrix elements still weakly depend on the depth of the 
optical potential via {B) but as J is small they can be 
approximated neglecting the i?-term and setting Jq = 1 
in the field expectation value a" = ri/[n~ i{Ac — NUq)]. 



u 



B 



(9) 




FIG. 2: (color online), (a) Numerically found average poten- 
tial depth Vcs = Vci -I- fiUo{a^ a) of the ground state of Q for 
N = 2 atoms in two wells, (b) Relative error in this potential 
using expansion @. The parameters are rj — 2k, Ki ~ ~4:Er 
and the scattering length is Us = OAEr. 



Lets now discuss some key physics. In the simplest 
case of a single particle in two wells the symmetric and 
antisymmetric superpositions of the atom in cither site 
are eigenstates with an energy difference 



AE = 2 



E + jlVci- hUoTf 



A' 



(10) 



strongly depending on the cavity parameters (see 
Fig. |3^). Hence detuning gives a handle to control the 
tunnel coupling and atom confinement (Fig. IJh). Note 
that the symmetric and antisymmetric eigenstate are as- 
sociated with different field amplitudes (lattice depths). 

Adding more atoms the interaction term comes to play 
and the ground state of the system is a superposition of 
different atomic configurations. Here the cavity param- 
eters influence the position and shape of the well known 
Mott-insulator superfluid transition |M Q - An impor- 
tant quantum feature appears for fields where the un- 
certainty in the photon number is not neglible. For an 
average photon number n generating a potential depth 
close to the phase transition point the photon numbers 
nil are than associated to different atomic phases, so 
that the ground state contains atomic states of different 
phases correlated with the corresponding photon num- 
ber. Even for a system being dominantly in the insulator 
phase, photon number fluctuations then allow the atoms 
to jump. This is shown in Fig.^Ji for 4 particles in 4 wells. 
Here we compare the site occupation probabilities as a 
function of scattering length for a purely classical and 
a quantum potential where the photon number uncer- 
tainty allows hopping even in the insulator regime. This 
behaviour can be enhanced or reduced through cavity 
mediated interaction as shown in Fig. ^p, where we plot 
the atom number fluctuations in one well as a function of 
flg for different atom-cavity detunings and compare it to 
the classical fleld case. Clearly the atom number fluctu- 
ations are enhanced on one side of the cavity resonance 
and suppressed on the other. 

The appearance of long range 4 particle interactions 
mediated by the i?^-term in Q can be seen by cornparing 
the density-density correlation functions {niUj ) M for 
nearest and next nearest neighbor sites. Depending on 
cavity parameters as shown in Fig. \^ each of the two 
correlations can be enhanced or reduced with respect to 
the classical potential case. 

Coupled atom-field dynamics in the semiclassical limit: 
As final point we turn to the classical field limit of the 
coupled Hamiltonian 10 for fixed N and large photon 
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FIG. 4: (color online), (a) Occupation probabilities pi for 
i atoms in one well as function of scattering length for the 
ground state of 4 atoms in 4 wells. The parameters are 
Uo — and Ac = —3.75k and Veft = hUo{a''a) — ~-4Er. 
For comparison the dashed lines correspond to an equivalent 
classical potential, i.e., (a^a) — and Ki ~ — 4£'j{. (b) Fluc- 
tuations of atom number in one well as a function of for 
different cavity detunings Ac = —5k (solid line), Ac — —3k 
(dashed line) and a classical field (dashed-dotted line). 
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FIG. 5: (color online), (a) Difi'erence of the density-density 
correlation functions: {runs} — (nin2). Parameters are as in 
Fig-Ela. (b) Dynamical evolution of atom number fluctuations 
in one well starting from an interaction free ground state and a 
sudden turn on of onsite interaction. A linear fit is depicted by 
the dashed line. Here we chose Uo — — k, Ac = —4.2k, Vd ~ 
0, as — SEr and ri such that Vcff = —4:Er. 



number. Here the field is approximately in a coherent 



state I a) and the system is assumed to evolve as a product 
state 1$) = 1-0) (8) \a:{t)). The Heisenberg equation (jZj) 
for the field then reduces to a classical equation for a 
containing expectation values of atomic operators: 



i (Ac - Uo{ij\JoN + JI3\Tp)^ - n a{t) + ry 



[E + j(v,i + hUo\a{t)\'] Bl^P) 



(11) 



a(t) is then inserted back into the atomic Hamilto- 
nian like a classical time dependent potential Vd 0, ^3 • 
Similar to the case of a time dependent Gross-Pitaevskii 
equation the corresponding Schrodinger equation can 
be solved simultaneously, where the matrix elements E, J 
have to be recalculated in each time step. 

Although the factorizing assumption is in general 
doubtful and one has to check the dynamical restriction 
to the lowest band, this procedure gives a first insight in 
the dynamical behavior of the model. As a generic ex- 
ample we show the time evolution of the uncertainty of 
the site occupation number starting with a 'superfluid' 
state at i = 0, when the onsite interaction is turned on. 
Fig. [SJd shows that in contrast to a fixed external poten- 
tial the dynamic cavity field leads to a damping of the 
fluctuations approaching a Mott insulator state. 

In summary we have shown that a dynamical quantum 
optical potential for ultracold atoms invokes a wealth of 
new physics. The effects arc pronounced in the limit of 
strong coupling and small photon numbers but long range 
interactions persist even in the bad cavity limit within a 
classical field approximation. The considered systems are 
in range of current experimental capabilities and should 
allow to control and study new quantum phases. 
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